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GLOBAL REGULARITY OF THE NAVIER-STOKES EQUATION 
ON THIN THREE-DIMENSIONAL DOMAINS 
WITH PERIODIC BOUNDARY CONDITIONS 

Stephen Montgomery-Smith 
Abstract 

This paper gives another version of results due to Raugel and Sell, and similar 
results due to Moise, Temam and Ziane, that state the following: the solution of the 
Navier-Stokes equation on a thin three-dimensional domain with periodic boundary 
conditions has global regularity, as long as there is some control on the size of the initial 
data and the forcing term, where the control is larger than that obtainable via "small 
data" estimates. The approach taken is to consider the three-dimensional equation as a 
perturbation of the equation when the vector field does not depend upon the coordinate 
in the thin direction. 

§1. Introduction 

The celebrated Navier-Stokes equation is concerned with the velocity vector u 
on a domain f2, describing the flow of an incompressible fluid. A famous unsolved 
problem is the following: if 17 is a nice enough 3-dimensional domain, and if the 
initial data is smooth, and the forcing term is uniformly smooth in time, then does 
it follow that the solution is smooth for all time? What is known is that a weak 
solution exists, although it is not known if that solution is unique. For details, we 
refer the reader to a number of books and papers, including [CF], [DG] and [T]. 

For the 2-dimensional problem, the solution is known. A heuristic approach to 
solving the 3-dimensional problem is as follows: if the solution becomes less smooth, 
then since we are dealing with an incompressible fluid, the complicated activity is 
going to get squashed into flat sheets, and one might hope that the solutions on the 
flat sheets can be somehow dealt with by the 2-dimensional case. Certainly, this 
"flattening" is observed in numerical and physical experiments. 

For this reason, it would seem that in order to get some handle on the real 
problem, that it might be important to understand what happens to the solution to 
the Navier-Stokes equation on thin domains, that is, domains of the form M x [0, e], 
where M is some 2-dimensional manifold, and e > is a small number. This is 
what Raugel and Sell did in a series of papers [RSI], [RS2], [RS3], as did Avrin [A], 
Temam and Ziane [TZ1], [TZ2], Moise, Temam and Ziane [MTZ], and Iftimie [II], 
[12]. 

In this paper, we consider the same situation considered by Raugel and Sell 
in [RS2], or by Moise, Temam and Ziane in [MTZ], or by Iftimie in [II], [12]. Let 
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Q e = [0, l\] x [0, 12] x [0, e], where l\ > I2 are some positive numbers, and e G (0, fa/£) 
is some small number. Let us consider vector fields u : O e — > R 3 satisfying the 
periodic boundary conditions 

«(0,y, z) = u(h,y,z), u(x,0, z) = u(x,l 2 , z), u(x,y,0) = u(x,y,e). (1) 

Given any appropriately smooth (where being in L2 is smooth enough) vector field 
u satisfying these boundary conditions, we may split it into its divergence-free part 
Lu and its gradient part. Thus L is the so called Leray projection. 

The Navier-Stokes equation considered in this paper is the equation for a func- 
tion u(t), t > 0, taking values in the space of 3-dimensional divergence- free vector 
fields on Q e satisfying the boundary conditions (1). The equation is 

d t u = vAu-L(u-Vu) + L(f), (2) 

where v is a positive constant (the viscosity), and / is a function of t taking values 
in the 3-dimensional vector fields satisfying the boundary conditions (the forcing 
term). For simplicity let us suppose / = L(f). 

In fact, to simplify our presentation, it will make sense to assume that 
In f ^ = 0> an d that J n u(0) dV = 0. In that case it is not hard to see that we have 
that J n u(t) dV = for all i > 0. This assumption does not really affect the solution 
very much, for suppose that we do not have this assumption. Given any function g on 

let g denote its mean value (hhc) -1 f n g dV . Let ($,t,Vt, Ct) = J u(0) + f(s) ds. 
Then replacing u(x, y, z, t) by u(x + £t, y + rj t , z + (t,t) — u gives us a solution to the 
Navier-Stokes equation in which / is replaced by / — /, and in which u = 0. 

It is known that in order to show global regularity of u that it is sufficient to 
show that u stays within the Sobolev space H 1 , that is, the space in which one 
derivative is in L 2 - Furthermore, once this is established, it also follows that the 
solution is unique. (See [CF] or [DG].) We will also include results concerning the 
Sobolev space H 2 , the space in which two derivatives are in 

Throughout this paper, we will use the letter c to denote a positive constant 
(typically larger than one), whose value will change with each occurrence. Only 
in certain places (such as in Lemmas 3 and 5) will we use subscripts on the c's to 
identify them, so as to avoid confusion. 

Theorem 1. Let u satisfy the Navier-Stokes equation (2) with periodic boundary 
conditions (1), and f Q udV = 0. Let 

U= \\u(0)\\ H1 , F = sup||/(t)|| 2 , M = max 



1^ 
h 

u(t) is in H 1 for all t > 0, with 



vl 

If M < c~ x —j — , then there exists a solution with the following properties. First, 



{ Z 3/2 

\u(t)\\ m < cmax I M, e _1/2 M 2 
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If t > c— , then 



v 




e -l/2 F 2 



Second, u{t) is in H 2 for almost all t>0, and indeed 




u(s)\\ 2 H 2 ds 



< oo 



for all < t < oo. 

Using a rescaling argument, it may be shown that is is sufficient to show Theo- 
rem 1 in the case that li, I2 and v are of order one, say, that these quantities all lie 
between 1/2 and 2. We will demonstrate this at the end of the proof of Theorem 1. 
So for the remainder of this discussion, let us focus on this case. Then in effect, 
M = max{[7, F}, and Theorem 1 gives global regularity in the case that M < c _1 . 

This result is not obtainable by the usual "small data" arguments, because 
these would only give regularity in the case that M < c~ l e 1 ' 2 . (The e 1 / 2 comes 
from the fact that we are calculating an L 2 norm on a domain whose total measure 
is of order e.) 

We will obtain Theorem 1 by considering it as a perturbation of the Navier- 
Stokes equation in which neither u nor / depend upon the third coordinate of £l e 
(the V coordinate). This approach was also taken by Iftimie [II], [12], and by Moise, 
Temam and Ziane [MTZ]. Note that the following result effectively does not depend 
upon e, in that if the result is obtained for one value of e, it is then automatically 
obtained for all others by a rescaling argument. 

Theorem 2. Let u satisfy the Navier-Stokes equation (2) with periodic bound- 
ary conditions (1 ), and J n u dV = 0, and neither u nor f depend upon the third 
coordinate of $7 e . Let 



Then there exists a solution with the following properties. First, u{t) is in H 1 for 
all t > 0, with 



17= ||u(0) || 




t 



2' 





Ift > c— , then 



v 




Second, u(t) is in H 2 for almost all t>0, and indeed 
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In order to compare Theorem 1 with the results in the literature, let us define 
the following projections. Let 

Pu(x, y, z ) = - I u(x, y, C) d(, 



and let Qu = u — Pu. As we stated above, we will discuss these results only in the 
case when Zi, I2 and v are of order one. Then the result of Raugel and Sell [RS2] 
gives global H 1 boundedness of the solution as long as 

||P^(0)|| H1 <e^ +5l (log(lA)) 52 

\\Qu(0)\\ H1 <e-^(l°g(lA)) 54 

sup||P/(t)|| 2 <6^(log(l/e)) 56 
t 

sup ||Q/(t)|| 2 (log(1/e)) * 
t 

where <5j (1 < i < 8) are positive numbers. 

The result of Moise, Temam and Ziane [MTZ] gives H 1 boundedness for t G 
[0, T(e)], where lim e ^ T(e) = 00, and also an integral condition on the H 2 norm, 
as long as 

\\Pu(0)\\ H1 <a(e)e L e+ s 
\\Qu(0)\\ H ,<a(e)e- k e+ s 

sup||P/(t)|| 2 <a(€)eH« 
t 

sup||Q/(t)|| 2 <a(e)€-H* 
t 

where 5 is a positive number, and lim e ^o a ( e ) = 0. 

Iftimie [II], [12] gets global existence results under conditions that use certain 
'anisotropic' Sobolev spaces. For example, one case of his Theorem 4.1 gives global 
existence under the condition that the forcing term / is zero, and 

||Qw|| H i/ 2 exp(ce -1 ||Pii||2) < c _1 . 

Even though his conclusions are slightly different, it is instructive to see how his 
hypothesis relates to that of this paper. Indeed, his condition is true if 

WPuW^^c^e^y^iiTJe) 
\\Qu\\ H1 <c- l e- l / 2+5 , 

where c$ depends upon 5 > 0. 

Before proceeding further, let us set our notation, and also state the "tools of 
the trade," that is, the standard results that many people in this area use. (See [CF] 



EJDE-1999/11 Navier-Stokes equation 5 

or [DG] — in particular [DG] considers the case of periodic boundary conditions, 
and many of the following calculations may be found there.) 

As well as the operators P and Q given above, let us define the operators R 
and S: R((ui,u 2 ,u 3 )) = (ui, u 2 ,0), and Su = u — Ru = (0, 0, u 3 ). Split u = v + w = 
r + s + w, where v = Pu, w = Qu, r = Rv, s = Sv. Since r and s do not depend 
upon z (where we label the coordinates of Q e by x, y and z), it is clear that r, s, v 
and w are all divergent free vector fields. 

As a notational device, I will write 



for the Sobolev space a derivatives in L p . (Thus D might represent the operator 
v 7 — A.) We have the Sobolev inequalities: if / is a mean zero function on [0, Zi] x 
[0, 12], 1 < p < q < 00, and a > 0, then 

11/11, <4D a f\\ p , 

where - + ^ = ^, and c depends upon p and q, as well as l\ and l 2 . Thus, if / is a 
mean zero function on Q e , then 

IIP/ll^ce-^Hzrp/i^. 

In fact, the only condition under which we will use this inequality is in the case 
p = 4 and q = 2, when a = 1/2. For this case, we will include an elementary proof 
in the Appendix. 

We have the interpolation inequality: if / is a mean zero function, ao,ai are 
real numbers, and < 9 < 1, then 

\\D^f\\ 2 <c\\D^f\\\- e \\D^f\\l 

where ag = (1 — 9)ao + 9a\. This inequality is easy to show using Parseval's 
identity and Holder's inequality. (See for example the proofs of Lemmas 4 or 6 for 
the statement of Parseval's identity.) 

We also have the Poincare inequalities: if / is a mean zero function on f2 e , with 
periodic boundary conditions, and a > 0, then 

Il/I|2<CP Q /Il2. 

where c depends upon a as well as l\ and / 2 . Again, this is easy to show using 
Parseval's identity. 

If u is a divergence-free vector field on the domain Q e with periodic boundary 
conditions, and if / and g are two other functions on £l e with periodic boundary 
conditions, sufficiently smooth so that the following integrals make sense, then by 
integration by parts we get 

/ f(u-Vg)dV = - f g(u-Vf)dV, 
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and so 

f(u-Vf)dV = 0. 



J 



If r is a two-dimensional, divergence-free vector field on the domain [0, li] x [0, l 2 ] 
with periodic boundary conditions that is sufficiently smooth, then we have the 
"enstrophy miracle:" 

/ Ar • (r • Vr) dA = 0. 

J[o,h]x[OM 

This is obtained as follows. First, integrating by parts, we see that the left hand 
side is equal to 

- / d x r ■ {8 x r • Vr) dA - [ d x r • (r • Vd x r) dA 

J[o,h]x[o,i 2 } ^[o,/i]x[o,; 2 ] 

— / d y r ■ (d y r ■ Vr) dA — d y r ■ (r • Vd y r) dA. 

J[o,h]x[o,h] ' J[0,h]xlO,h] 

(Here, as in the rest of the paper, d x , d y and d z represent partial differentiation with 
respect to x, y and z respectively, that is, the first, second and third coordinates 
respectively.) We see that the second and fourth terms are zero. Expanding and 
collecting the first and third terms, and remembering that r is divergence-free, we 
see that they also total to zero. 



Navier-Stokes for Flows Independent of z 

Let us start with the proof of Theorem 2, the case when w = Qu = 0, that is, when 
u = Pu = v. We will prove Theorem 2 in the case that the quantities Zi, 1% and v 
all lie between 1/2 and 2. The general result may be obtained as shown at the end 
of the proof of Theorem 1. 

First we need to find a solution to the Navier-Stokes equation. This is done 
using so called Galerkin solutions. Let S n denote the projection that takes a function 
/ on i7 e onto the nth partial Fourier series. (Quite how we index this sum is not 
important, as long as S n converges formally to the identity.) Then we consider the 
solution u n to the problem 



d t u n = uAu n - S n L(u n • Vu n ) + S n L(f), 



with u n (0) = S n u(Q) for which u = S n u. It is a well known argument to show 
that this equation, essentially an ODE, has a global solution, and that the solutions 
u n converge weakly to some function u. This is the so called weak solution to the 
Navier-Stokes equation. In that case, for any appropriate norm ||-||, we will have 
that ||n|| < liminfn^oo ||n n ||. 

Thus, in order to prove our theorem, it is sufficient to prove it for the Galerkin 
solutions. This is what we shall do. However, carrying the symbol S n throughout 
the proof could be a little cumbersome, and so for this reason, we will replace L by 
S n L, and suppose that both / and u(0) lie in the range of S n . (We will also suppose 
that / lies in the range of L.) 
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d t v = uAv — L(r ■ Vv) + /, 



since v ■ Vv = r • Vv, because d z v = 0. If we apply the R and S operator to this, we 
get the following pair of equations: 



The first equation is merely the 2-dimensional Navier-Stokes for the flow r. The 
second equation essentially says that the 1-dimensional quantity s is being pushed 
around by the 2-dimensional flow r (and indeed in the second equation, the operator 
L acts as the identity). 

Let us write 

4> = \\Dr\\ 2 , ip = \\Ds\\ 2 , 4> = \\D 2 r\\ 2 , ip = \\D 2 s\\ 2 , 9 = \\u\\ 2 . 

Poincare's inequality tells us immediately that 4> < afi, ip < cip, and 9 2 < c(4> 2 + ip 2 ). 

The process for comprehending cf) and (f> is well known. Start with the first 
equation in (3), dot product both sides with — Ar, integrate over O e , use the self- 
adjointness of L, apply some integration by parts, use the Cauchy- Schwartz inequal- 
ity, and use the "enstrophy miracle," to get 



Use the inequality ab < a 2 + b 2 to get that ||L> 2 r|| 2 F < (v /2)\\D 2 r\\ 2 2 + {2/u)F 2 . 
Thus we have that 



This differential inequality is easy to solve, but before we do so, let us first understand 
ifr and ip. Take the second equation from (3), dot product both sides with —As, 
integrate over and work as before. But in this case, the "enstrophy miracle" 
does not work — there is a term: 



d t r = vAr — L(r • Vr) + Rf 
d t s = uAs - L(r • Vs) + Sf. 



(3) 




d t <p 2 < -c" 1 ^ 2 + cF 2 . 




As • (r • Vs) dV. 



To get a grip on this term, see that it splits into 




plus another term with d 2 in place of d 2 . 



Integrate by parts to get 



f d x s ■ (d x r ■ Vs) dV - f d x s ■ (r • Vd x s) dV. 
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The second term is zero. For the first term, we may use Holder's inequality and the 
Sobolev inequality to bound it by: 

||i>|| 2 ||r>s||4 < ce- 1,2 \\Dr\\ 2 \\D 3/2 s\^ 2 

< ce _1/2 ||Z)r|| 2 ||L>,s|| 2 ||L> 2 ,s|| 2 

< ^||L> 2 s|| 2 + ce- 1 || J Dr||2||Ds|| 2 , 

where in the last step we use the inequality ab < a 2 + b 2 . 

Putting this all together, we get a differential inequality: 
d t 4> 2 < -c"V + ce" + cF 2 . 

We will also require a differential equation for 9: take the Navier-Stokes equa- 
tion, dot product both sides with u, integrate over O e , and do the usual stuff, to 
get 

-<9t||u|| 2 < — ^ll-D^I^ + l|£MI 2 -F- 



Since ||D«|| 2 F < (is/2)\\Du\\ 2 + (2/v)F 2 , we get 

-^<9 t ||d| 2 < -J||Du||o + -F 2 , 

that is 

dt e^<- c -\^ + ^) + C F 2 . 



Then Theorem 2 will be established once we have obtained the following result. 

Lemma 3. Let U, F and T be positive numbers, where T may be infinity. Let <f>, 
ip, <f>, tp, be positive differentiable functions of t. Suppose that for some positive 
constants Ci (1 < i < 10) we have 

0(0) < U (4.1) 

^(0) < U (4.2) 

2 < c 1 (<p 2 + ij 2 ) (4.3) 

<t> < c 2 4> (4.4) 

%j> < c 3 4> (4.5) 

dttf < -cl 1 ^ 2 + c 5 F 2 (4.6) 

dti? < -c^ip 2 + c 7 e~ W + c 8 F 2 (4.7) 

d t 6 2 < -eg 1 {4? + V 2 ) + c w F 2 (4.8) 

for < t < T. Let M = max{U, F}. Then there exist positive constants Ci 
(11 < i < 17), depending only upon a (1 < i < 10) such that we have the 
inequalities 

2 <cn(F 2 + (U 2 - F 2 )e- C ^ 1 ) (4.9) 

<P 2 < c 13 (F 2 + (U 2 - F^e' ™*) (4.10) 
$ < c 15 max{e~ 1/2 M 2 ,M}, (4.11) 

and if t > C\q then 
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(4.12) 



for < t < T. Furthermore, we have that 

t 

(0(s) 2 + V(s) 2 ) ds < oo (4.13) 
for < t < T. 

Proof: Inequalities (4.9) and (4.10) are easy to obtain from combining (4.1) and (4.4) 
with (4.6), and (4.1), (4.2) and (4.3) with (4.8), by using Gronwall's inequality. 

Let us obtain (4.11). From (4.10), we see that there is a positive number c\s 
such that <p < C18-M. Combining this with inequalities (4.5), (4.7) and (4.8), we see 
that for some positive constants cig and c 2 q that 

d t ^ 2 + c 19 V> 2 < c 20 e- 1 M 2 (F 2 - d t 6 2 ) + c 20 F 2 . 

Multiply both sides by e Cl9 * and integrate from to t to get that 

e ci9ty2 _ ^( )2 < C2() f e -l M 2 e ci9 S ( F 2 _ ds Q( s ^ + F 2 e c 19 s ^ 

Jo 

which, by evaluating the integrals, and integrating by parts, is less than or equal to 

c 2 oci 9 1 e Cl9t (e- 1 M 2 F 2 +F 2 )+C2 e~ 1 M 2 (e(0) 2 -e Cl9t e 2 )+c 20 c 19 e- 1 M 2 [ e Cl9S 9(s) 2 ds. 

Jo 

Now, from (4.9), we see that there is a positive constant c 2 i such that 9 < c 2 \M. 
Hence 

e Cl9 V 2 - ^(0) 2 < c 20 ci 9 1 e c ^ t (e- 1 M 2 F 2 + F 2 ) + c 20 c 2 l e- 1 M A + c 20 c 19 e~ l M 2 'c 2 21 f e Cl9S ds 

Jo 

< c 20 c^ 9 1 e Cl9t (e- 1 M 2 F 2 + F 2 ) + c 20 c 2 2l e- 1 M A + c^ 1 M 2 c 2 2l e c ^ . 

Hence 

V> 2 < e- Cl9t C/ + c 20 c^ 9 1 (e- 1 M 2 F 2 + F 2 ) + c^c^e"^-^ 9 ^ 4 + c^c^e^M 2 , 

and from here it is easy to obtain (4.11). 

To obtain (4.12) is similar. We see that there are positive numbers c 22 , c 2 s and 
C24 such that if t > c 22 then <f> < c 2 ^F and 9 < c 2 ±F. Apply the above argument, 
except integrate from c 22 to t instead of from to t. 

Finally, (4.13) may be obtained by integrating (4.6) and (4.7), and using (4.9), 
(4.10) and (4.11). Q.E.D. 
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Proof of Theorem 1 

The proof of Theorem 1 will follow the same lines as the proof of Theorem 2, with 
some additional work for dealing with the w = Qu part. Let us start by assuming 
that the quantities Zi, h and v all lie between 1/2 and 2. 

We need a couple of Poincare/Sobolev type inequalities on f2 e . 
Lemma 4. Let w = Qu, then 

IML ^ Cel/2 \\ D2w \\ 2 

IML ^ cel ^ 4 ll-^Hl2- 

Proof: Let w denote the Fourier coefficients of w: 

w(m, n,p) = (Zi^e) -1 / w(x,y, z) exp(—2iri(mx/li + ny/h + pz/e) dx dy dz, 

where m, n and p are integers. Then the function w can be reconstructed using the 
Fourier series 

w(x,y,z)= w(m,n,p)exp(27ri(mx/li + ny/l 2 + pz/e). 

m,n,p 

We recall Parseval's identity: 

\\w\\ 2 = ZiZ 2 e ^2 \w(m,n,p)\ , 

m,n,p 

and the Hausdorff- Young inequality: if 2 < p < oo, and p' = p/{p — 1), then 

i/p> 

w\ 



\m,n,p J 



Since u> = Qu, we have that u>(m, = if p = 0. 
Now for any real number a, we have that 

D^w(m, n,p) = {-2iti) a {m 2 /l 2 + n 2 /l 2 + p 2 /e 2 )" /2 ^K n,p), 
and thus by Parseval's identity 

\\D a w\\ 2 = (2n) a (l 1 l 2 e) 1 / 2 ( Yl ( m V*? + ^ + P 2 /f 2 )>(™> «,P)| 2 ) 



\m,n,p 



Let us start with showing the first inequality. Apply Cauchy-Schwartz to get 

< /p#0 ( m 2 A 2 +n 2 A 2 +p2/e2) -2 > J 7 x 

\m,n,p / 

( (m 2 /l 2 +n 2 /l 2 +p 2 /e 2 ) 2 \w(m,n,p)\ 



v 1/2 
2 I 



\m,n,p 
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By approximating sums by integrals, and using other elementary inequalities, we see 
that 

V I p + Q (m 2 /ll +n 2 /l 2 2 +p 2 /e 2 )- 2 < c V I p¥=Q . 1 . . . 

^ p ^ y 11 ' z ' ' ^ p ^ max{m 4 ,n 4 ,p 4 /e 4 | 



m,n,p m,n,p 



max{n 3 , p 3 /e 3 } 



J 9/9 
p c < 



<ce 2 . 



Hence we obtain the first inequality. 

The second inequality has a similar proof: start by using Holder's inequality to 



get 



\w\\ 4 <(hl 2 e) 1/4 (j2 M^,n, P )\ i/3 ) 

\m,n,p / 

< (/i/ 2e ) 1/4 ( £ IpM^/li + nVll+p'/e 2 )- 2 ] 

\m,n,p / 

( E ( m2 A? + « 2 A 2 2 +p 2 A 2 )K-»,«,p)I 



1/4 

X 



1/2 

2 I 



\m,n,p 



and proceed as with the proof of the first inequality. Q.E.D. 

Proof of Theorem 1: As in the proof of Theorem 2, we argue that we work with the 
Galerkin approximations. We will obtain differential inequalities for the following 
quantities: 

4> = \J \\Dr\\ 2 2 + |^Du^, ip = \J i^Dsjl^ + H-Dio^, 

X = \\D 2 w\\ 2 , 9=\\u\\ 2 . 

Poincare's inequality tells us immediately that <p < cifi, ip < dp, and 6 2 < c((p 2 + ip 2 ). 

Let us start with the Navier-Stokes equation, and apply the operator P. Note 
that if / and g are functions on Q e , then P((Pf)(Qg)) = 0. Thus, we obtain 



d t v = vAv - LP(v • Vu) - LP(w • Vw) +Pf = vkv - LP(r • Vu) - LP(w • Viu) + P/. 
Now apply i? to both sides: 

d t r = uAr - LP(r ■ Vr) - LP(w ■ VRw) + PRf. 
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Finally, take the dot product of both sides with — Ar, and integrate over O e , and do 
all the usual stuff. A lot of the terms work in exactly the same way that they did 
in the previous section. The only term that we did not deal with is the following: 



This splits into a term: 



/ Ar ■ (w ■ VRw) dV. 
Jn e 

[ d 2 x r ■ (w ■ VRw) dV, 



and a similar one with d 2 in place of d x . The bounds for the second term will be as 
for the first, so let us only deal with the first. Integrate by parts to get 

- / d x r ■ (d x w ■ VRw) dV - [ d x r ■ (w ■ Vd x Rw) dV. 
The first term is bounded by 



and the second by 



1 1 2 

c||-Dt>|| 2 ||Z)ti;|| 4 , 



c||Dt;|| 2 ||D HUMloo; 



Combining all this with Lemma 4, we get 

dt\\Dr\\ 2 2 < -c^WD^Wl + ce 1/2 ||£>u|| 2 ||D 2 «;||2 + cF 2 . (5) 

The work for s is practically identical, and we get 

dt\\Ds\\ 2 2 < -c- 1 \\D 2 s\\ 2 2 + ce- 1 || J Dr||2|| J Ds||2 + ce 1/2 \\Dv\\ 2 \\D 2 w\\l + cF 2 . (6) 

We also need to establish an equation for w. Take the Navier-Stokes equation 
and apply Q. Note that if / and g are two functions on fi e , then Q((Pf)(Pg)) = 0. 
Thus we get 

d t w = i/Ai« - LQ(w ■ Vv) - LQ(v ■ Vw) - LQ(w ■ Vw) + Qf. 

Take the dot product with — Aw, and integrate over f2 £ , doing all the stuff as before. 
Let us see what happens to the non-linear terms, only bothering with the d 2 w part 
of Aw, knowing that the other parts will give the same estimates. 

First we get 

/ d 2 x w ■ (w-Vv)dV = - [ d x w ■ (d x w ■ Vv) dV - [ d x w ■ (w ■ Vd x v) dV. 
Jn e Jn e Jn e 

The first term is bounded in absolute value by 

||D«;||4l|£)t;|| 2 < ce 1 / 2 \\D 2 w\\ 2 2 \\Dv\\ 2 . 
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The second term is equal to 



/ d x v-{ 



w ■ Vd x w) dV, 



which is bounded in absolute value by 



IMUI^HUI-^I^ ^ cel/2 \\ D2w \\l\\ Dv \\2- 

Next, we have 

/ d 2 w ■ (v ■ Vw) dV = — I d x w ■ (d x v ■ Vw) dV — / d x w ■ (v ■ Vd x w) dV. 
Jn e Jn e Jn e 

The first term is bounded in absolute value by 

||D«;||4||£)v||2 < ce 1/2 \\D 2 w\\l\\Dv\\ 2 , 

and the second term is zero. 
Finally we have 

/ d 2 w ■ (w ■ Vw) dV = — I d x w ■ (d x w ■ Vw) dV — / d x w ■ (w ■ Vd x w) dV. 
Jn e Jn e ifi £ 

The first term is bounded in absolute value by 

\\Dw\\ 2 \\Dw\\ n < ce 1/2 ||I} 2 w|| 2 ||LML, 



and the second term is zero. 

So, doing all the same stuff as above, we get 

0t\\Dw\\ 2 2 < -\\D 2 w\\ 2 2 + ce 1/2 \\Dv\\ 2 \\D 2 w\\ 2 2 + ce 1/2 \\Dw\\ 2 \\D 2 w\\ 2 2 + cF 2 . (7) 

If we add equations (5) and (7), and also (6) and (7), (and apply liberally the 
inequalities \J a 2 + b 2 < a + b < y/2 v 7 a 2 + b 2 for positive a and b) we get the two 
differential inequalities: 

M 2 < -c- l 4> 2 - c- V + ce x '\<f> + ^) X 2 + cF 2 

d t ^ 2 < -c- 1 ^ 2 - c- V + ce-VV 2 + ce 1/2 (0 + VOx' + cF 2 . 

In addition, arguing as in the previous section, we get the differential inequality 

^ 2 <-c-V+V> 2 )+cF 2 



Thus the theorem will be established when we have proved the following lemma. 
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Lemma 5. Let U and F be positive numbers. Let 4>, ip, <p, ip, 6 be positive 
differentiable functions oft. Suppose that for some positive constants a (1 < i < 10, 
18 < i < 19) we have 

0(0) < U (8.1) 

iP(0) < U (8.2) 

^ 2 <ci(^ 2 +V 2 ) (8.3) 

<p> < c 2 <t> (8.4) 

ip < c 3 ip (8.5) 

dttf < (-cfg 1 + c^e- 1 ' 2 ^ + if))) X - cjV 2 + c 5 F 2 (8.6) 

<9 t V 2 < (-crs + ci 9 e" 1/2 (0 + V0)% ~ ce 1 ^ 2 + c 7 e" VV + c s F 2 (8.7) 

d t 9 2 < -c^\4> 2 + ip 2 ) + c 10 F 2 (8.8) 

for < t < oo. Let M = max{U, F}. Then there exist positive constants Ci 
(11 < i < 17, i = 20), depending only upon Ci (1 < i < 10) such that if M < c^ 1 , 
then we have the inequalities 

2 <cn(F 2 + {U 2 - F 2 )e- C ^*) (8.9) 

2 <c 13 (F 2 + (U 2 -F 2 )e- c ^ t ) (8.10) 

ip < c 15 max{e~ 1/2 M 2 ,M}, (8.11) 

and if t > cie then 

ip < c 17 max{e~ 1/2 F 2 , F}. (8.12) 
for < t < oo . Furthermore, we have that 

f {4>(s) 2 + if)(s) 2 ) ds < oo (8.13) 
Jo 

for < t < oo. 

Proof: Let (ci)n<i<i7 depend upon (ci)i<i<io as in Lemma 3. Let us suppose that 
U,F < C20 , where C20 will be chosen momentarily. Let 

T = inf{t > : Ci 9 e 1 / 2 (0 + if)) > c^ 1 }. 

Suppose for a contradiction that T < 00. But then for t G [0,T], the quantities <p>, 
■ip, (ft, 4> and 9 satisfy the hypothesis of Lemma 3. But then by the conclusion of 
Lemma 3, we know that for some constant C21 > that 

+ if)< c 21 e- 1 / 2 cz 2 . 

Setting C20 small enough, we see then that for t G [0, T] that we have 

Cl9e i/2 (0 + v ,) < c -i/2, 

and thus that there is a neighborhood of T such that 

c 19 e 1 / 2 (0 + ^) < Cl8 \ 
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contradicting the definition of T. 

Thus T = oo, and thus the functions 4>, if), (ft, ip and 6 satisfy the hypothesis of 
Lemma 3, and the result follows. 



Now let us relax the restriction that h, l 2 and v all lie between 1/2 and 2. Let 

nlo 

if 



n be the integer part of -i, and define new vector fields u and / on [0, 1] x 



according to the formulae 



l\ ( I 2 \ 

u(x 1 ,x 2 ,x 3 , t) = — u( lixi, {hx 2 mod Z 2 ), hx 3 , -ji j , 

\\ ( I 2 \ 

f(x 1 ,x 2 ,x 3 ,t) = — f I l lXl , (hx 2 mod / 2 ), /1X3, ^-t J . 

Then it may be easily verified that these satisfy the equation 

8 t u = Au - L(u ■ Vu) + L(f), 

that is, one may apply the version of Theorem 1 that we already have to the func- 
tions u and /. Obtaining the more general version of Theorem 1 is then merely a 
question of interpreting what it says about u and /, taking into account the following 
identities: 



v 2 



2 ' 
n 



1/2J3/2II-MI2' 



Ffi — TTo" ^ ffi- 



Q.E.D. 
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Appendix: the Sobolev Inequality 

The following result is essentially part of the literature. For example, in [S] this 
result is found for functions on Euclidean space. However, for our special case, we 
are able to provide an elementary proof (a proof motivated by Littlewood-Paley 
theory) . 

Lemma 6. Let f be a function on [0,l\] x [0,l 2 ] satisfying periodic boundary con- 
ditions, that is mean zero. Then there is a positive constant c, depending only upon 
l\ and l 2 , such that 

\\fh<c\\D^f\\ 2 . 

Proof: For each r = (r 1 ,r 2 ) a pair of integers, write \r\ = \Jr\jl\ + r 2 /l 2 . Define 
the Fourier coefficients of /: 

rli rh 

fr = (hh)~ 1 I I f(x,y)exp(-2m(r 1 x/l 1 +r 2 y/l 2 ))dydx. 
Jo Jo 

The original function can be reconstructed using the Fourier series 
f{x,y) = ^/ r exp(2vri(r 1 x// 1 + r 2 y/l 2 )), 

r 

and we have Parseval's identity 

\\f\\ 2 2 = (hh)^\fr\ 2 - 

r 

We see that 

D^f r = >/=fcri(r?/l? + r 2 2 // 2 2 ) 1/4 /„ 
and so by Parseval's equality we see that 

\\D^f\\l = 2vr(M 2 ) £(r?/*? + rl/H)^\f r \ 2 . 

r 

For m a non-negative integer, set 

A m = J y ] \fr 

\2 m <\r\<2 m + 1 

Notice then that 

00 

£ < c\\D^f\\ 2 2 . 

m=0 

Now, 

\\ft = I' I 2 f(x,y) 2 7&Wdydx, 
Jo Jo 
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and expanding this using the Fourier series, we obtain that 



(4) 



4 — (hh) E frWfrWfrWfrW 

r (i)_|_ r (2)_ r (3)_ r (4) =0 

< 24(Zi/ 2 ) E |/ r (i)/r( 2 )/r( 3 )/r 

r (l) +r (2)_ r (3)_ r (4) =Q 

which in turn is bounded above by 

24(M 2 ) £ E E l^ 2) l 

0<m 1 <m 2 <m 3 2 m i < rW | <2 m i + 1 2 m 2 <|r< 2 ) |<2 m 2 + 1 

^] /| r (l) +r (2)_ r (3)|>| r (3)| |/ r (3)/ r (l) +r .(2)_ T .(3) |. 

2 m 3<|r< 3 )|<2 m 3 + 1 

In bounding this quantity, let us start by looking at the inner sum: 

/| r (l) +r (2)_ r (3)|>| r (3)| |/ r (3)/ r (l)_)_ r (2)_ T .(3) |. 

2 m 3<|r( 3 )|<2 m 3 + 1 

Using the Cauchy- Schwartz formula, this can be bounded above by 

( E l/r< 3 )| 2 ) ^ ( E ^| r (l)+ r (2)_ r (3)|>| r (3)| |/ r (l) +T ,(2)_ r (3) | 2 ) ^ . 

2 m 3 < |r( 3 ) | <2 m 3 + ! 2 m 3 < |r( 3 ) | <2 m 3 + x 

The first term in this product is A m3 . As for the second term, since mi < m,2 < m 3 , 
it follows that \r^ + r^ 2 ^ — r^\ < 3 • 2™ 3 , and hence the second term is bounded 
by (^4^3 + ^4^3+1 + ^4 2 n3 +2) 1 ^ 2 - Thus, the above quantity can be bounded above by 

Furthermore 

E \M<( E 1) ( E \fA <c2 m A 



m j 



2™<|r|<2 m + 1 \2"*<|r|<2 m + 1 / \2 m < |r | <2 m + 1 



since the number of points r such that 2 m < \r\ < 2 m+1 is 2 2m to within a constant 
factor. 

Thus 

OO m 3 7712 

H/H4 < C E (^i 3 +^m 3 + l + J 4m 3 4-2) E 2 m2 Am 2 E "^Arn- 
m 3 =0 m 2 =0 m 1 =0 

Now, applying Cauchy-Schwartz, we see that the inner sum obeys the inequalities 

TTL2 I TU2 \ V2 / m2 \ 1/2 

E 2 " ll ^ 1 < E H E 2mM ^ ^ c2^/ 2 || J D 1 /2 / || 2 . 

mi=0 \mi=0 / \mi=0 / 
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Thus 

oo m 3 

ll/lll < 4D^f\\ 2 £ « 3 + A 2 m3+1 + < 3+2 ) £ 2 3 -/ 3 A m2 . 

m3=0 1712=0 

Applying Cauchy- Schwartz once more in a similar fashion, we get that 

Y, 2 3m ^ 3 A m2 <c2 m ^\\D 1 / 2 f\\ 2 , 

m 2 =0 

and so 

oo 

< 4D 1/2 f\\l E 2 m3 (^ 3 + Al 3+1 + Al 3+2 ), 

m 3 =0 

from which we see that H/H4 < cWD 1 / 2 /]^ as required. 

Q.E.D. 

Acknowledgments 

This paper owes a lot to other people who have helped me over the years. I would like 
to thank to Joel Avrin and John Gibbon who explained the Navier-Stokes equation 
to me three years ago. I would like to thank George Sell, Mohammed Ziane and 
Benoit Desjardins for improvements and corrections to this paper. 

References 

A Avrin, Joel D. Large-eigenvalue global existence and regularity results for the 
Navier-Stokes equation. J. Differential Equations 127 (1996), no. 2, 365-390. 

CF Constantin, Peter; Foia§, Ciprian, Navier-Stokes equations. Chicago Lectures in 
Mathematics. University of Chicago Press, Chicago, IL, 1988. 

DG Doering, Charles R.; Gibbon, J. D. Applied analysis of the Navier-Stokes equa- 
tions. Cambridge Texts in Applied Mathematics. Cambridge University Press, 
Cambridge, 1995. 

11 Iftimie, Drago§ Les equations de Navier-Stokes 3D vues comme une perturbation 
des equations de Navier-Stokes 2D. (French) [The 3D Navier-Stokes equations 
seen as a perturbation of the 2D Navier-Stokes equations] C. R. Acad. Sci. Paris 
Ser. I Math. 324 (1997), no. 3, 271-274. 

12 Iftimie, Drago§ The 3D Navier-Stokes equations seen as a perturbation of the 
2D Navier-Stokes equations, preprint, available at 

http : //www. maths .univ-rennesl . f r/~ if timie/publi cat ions .html 

MTZ Moise, I.; Temam, R.; Ziane, M. Asymptotic analysis of the Navier-Stokes equa- 
tions in thin domains. Dedicated to Olga Ladyzhenskaya. Topol. Methods Non- 
linear Anal. 10 (1997), no. 2, 249-282. 

RSI Raugel, Genevieve; Sell, George R. Navier-Stokes equations on thin 3D domains. 
I. Global attractors and global regularity of solutions. J. Amer. Math. Soc. 6 
(1993), no. 3, 503-568. 



EJDE-1999/11 



Navier-Stokes equation 



19 



RS2 Raugel, G.; Sell, G. R. Navier-Stokes equations on thin 3D domains. II. Global 
regularity of spatially periodic solutions. Nonlinear partial differential equations 
and their applications. College de France Seminar, Vol. XI (Paris, 1989-1991), 
205-247, Pitman Res. Notes Math. Ser., 299, Longman Sci. Tech., Harlow, 1994. 

RS3 Raugel, Genevieve; Sell, George R. Navier-Stokes equations in thin 3D domains. 
III. Existence of a global attractor. Turbulence in fluid Rows, 137-163, IMA Vol. 
Math. Appl, 55, Springer, New York, 1993. 

S Stein, Elias M. Singular integrals and differentiability properties of functions. 
Princeton Mathematical Series, No. 30 Princeton University Press, Princeton, 
N.J. 1970. 

T Temam, Roger Infinite-dimensional dynamical systems in mechanics and physics. 
Second edition. Applied Mathematical Sciences, 68. Springer-Verlag, New York, 
1997. 

TZ1 Temam, R.; Ziane, M. Navier-Stokes equations in three-dimensional thin do- 
mains with various boundary conditions. Adv. Differential Equations 1 (1996), 
no. 4, 499-546. 

TZ2 Temam, R.; Ziane, M. Navier-Stokes equations in thin spherical domains. Op- 
timization methods in partial differential equations (South Hadley, MA, 1996), 
281-314, Contemp. Math., 209, Amer. Math. Soc, Providence, RI, 1997. 

Stephen Montgomery-Smith 
Math. Dept., University of Missouri 
Columbia, MO 65211, U.S.A. 
email: stephen@math.missouri.edu 
http: / /math, missouri.edu/~stcphen 



20 



Stephen Montgomery-Smith 



EJDE-1999/11 



Addendum: July 8 1999. 

On page 3, in both Theorems 1 and 2, the sentence that begins "Ift > c — " should 
read 

Furthermore 



Similarly, for Lemma's 3 and 5 (pages 8/9 and 14 respectively), the phrase that 
begins "and if t > c w then" and ends "for < t < T" (respectively "for < t < oo" ) 
should be replaced by 



For Lemma 5 the corresponding equation number is (8.12). 

The second to last paragraph of the proof of Lemma 3 (page 9) should be 
changed to 

To obtain (4.12) is similar. We see that given e > there are positive 
numbers r, C23 and C24 (where r depends upon e) such that if t > r then 
4> < C23-F + e and 6 < C24-F + e. Apply the above argument, except integrate 
from r to t instead of from to t. 




and 



limsup^ < C17 max{e 



1 / 2 F 2 , F}. 



(4.12) 



t— s-oo 



